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Abstract 

When Xt *C (It, where Xt is the de Broghe wavelength and (It, the distance 
of closest approach of thermal electrons, a classical analysis of the energy of 
a plasma can be made. In all the classical analysis made until now, it was 
^ . assumed that the frequency of the fluctuations uj <^ T (ks = h = 1). Us- 

o. 

0^ ■ ing the fluctuation- dissipation theorem, we evaluate the energy of a plasma, 

allowing the frequency of the fluctuations to be arbitrary. We find that the 

Oh; 

Q . energy density is appreciably larger than previously thought for many inter- 

c/3 . esting plasmas, such as the plasma of the Universe before the recombination 

a 

>• . era. 
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I. INTRODUCTION 

There have been many classical calculations of the energy of a plasma |l|-^. They are 
based on perturbation theory of an ideal gas, in terms of the plasma parameter g (which 
usually is a small value). The treatment, to the first order in g, is the Debye-Hiickel theory. 
However, in the calculations that have been made it is assumed that uj <^ T [ks = h = 1). 
This is a very strong assumption. For example, in our previous analysis @,§], we showed 
that only by not assuming u; -C T, is the blackbody spectrum obtained. 

We evaluate the energy of a plasma, studying the electromagnetic fluctuations in a plasma 
without assuming that u; ^ T. A plasma in thermal equilibrium sustains fluctuations of 
the magnetic and electric fields. The electromagnetic fluctuations are described by the 
fluctuation-dissipation theorem [^. 

The evaluation of the electromagnetic fluctuations in a plasma has been made in nu- 
merous studies 0. Recently, Cable and Tajima (see also 0) studied the magnetic field 
fluctuations in a cold plasma description with a constant collision frequency as well as for a 
warm, gaseous plasma, described by kinetic theory. 

Using a model that extends the work of Cable and Tajima [§, we study an electron- 
proton plasma of temperature 10^ — 10^ K with densities 10^^ — 10^^ cm'^. The condition 
for a classical analysis is that Ay < d^, where A^ is the de Broglie wavelength for a thermal 
electron and dr = e^/T, the distance of closest approach. This condition is satisfied for 
T < 3 X 10^ K and for the plasmas studied. 

In section |l| we recall the expressions for the electromagnetic fluctuations in a plasma, 
and in section |T|, the electromagnetic energy is computed. Finally, we discuss our results 



in section IV 



II. ELECTROMAGNETIC FLUCTUATIONS 



The spectra of the electromagnetic fluctuations in an isotropic plasma are given by 






■i-n e^i-" — I \ El 



---ft)" 






{h = ks = 1), where El and ^t are the longitudinal and transverse dielectric permittivities 
of the plasma. The first and second terms of Eq. (|1|) are the longitudinal and transverse 
electric field fluctuations, respectively. 

By using the fluctuation-dissipation theorem, we can estimate the energy in the electro- 
magnetic fluctuations for all frequencies and wave numbers. The calculation includes not 
only the energy of the fluctuations in the well defined modes of the plasma, such as plasmons 
in the longitudinal component and photons in the transverse component, but also the energy 
in fluctuations that do not propagate. 

For the description of the plasma, we use the model described in detail in Opher and 
Opher [Q^ . The description includes thermal and collisional effects. It uses the equation of 
Vlasov in first order, with the BGK (Bhatnagar-Gross-Krook) collision term that is a model 
equation of the Boltzmann collision term |]I0| . We used the BGK collision term as a rough 



guide for the inclusion of collisions in a plasma. 

From this description, the dielectric permittivities for an isotropic plasma are easily 
obtained: 
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where a is the label for the species of particles, v^ the thermal velocity for the species and 
Z{z), the Fried and Conte function. 



III. ELECTROMAGNETIC ENERGY 

In order to estimate the electromagnetic energy, we use the dielectric permittivities, given 
by Eqs. @ and (^), and calculate the magnetic and the electric field spectra from Eqs. (|l]) 
and (^. Integrating the spectra in wave number and frequency (and dividing by {2tt) ), we 
obtain the energy densities of the magnetic field pB and of the transverse and longitudinal 
electric fields pej- and pl. 

Usually, when estimating the energy stored in the electromagnetic fiuctuations from Eqs. 
(0) and (0), it is assumed that uj <^T [ks = h = 1). With this assumption, the Kramers- 
Kronig relations can then be used, and a simple expression for the energy is obtained [|l],^ . 
However, the assumption that u; ^ T is very restrictive. For example, a large part of 
the fiuctuations which create the blackbody electromagnetic spectrum has uj > T 0,§]. It 
is therefore necessary to perform the integration of the spectra over frequency and wave 
number without using this assumption. 

Our model uses kinetic theory with a collision term that describes the binary collisions 
in the plasma. A cutoff has to be taken since, for very small distances, the energy of the 
Coulomb interaction excedes the kinetic energy. This occurs for distances Tmin ~ e^/T, 
which defines our maximum wave number, k^ax- 

A large kmax is needed in order to reproduce the blackbody spectrum. In this study, we 
used a k^ax equal to the inverse of the distance of closest approach, which we previously 
found is able to do this |^,^ . Any smaller k^ax was unable to reproduce the entire blackbody 
spectrum. 

In the usual classical calculations, the correction to the energy due to correlations between 
the particles, is made through the correlation energy. To the first order in the plasma 
parameter g, the correlation energy depends on the two-particle correlation function S{k), 

Ec = ^j dkk^kS{k) - £^ / dkk^k , (5) 

where the second term is the energy of the particles due to their own fields. S{k) can be 
estimated by the fiuctuation-dissipation theorem or by the BBGKY hierarchy equations [§. 



Generally, it is assumed that uj <^T (so the Kramers- Kronig relation can be used) and S{k) 
is obtained as 



S(*) ^ ^ . (6) 

where ko is the inverse of the Debye length. 

With this, the energy density of a plasma to first order in g is given as 

where n is the number density of the particles. Thus, the correlation energy, to the first 
order in g is 

We define the energy of a plasma as 

f/=^nT(l + A). (9) 

With this definition, A = Aq = —g/12TT, for the previous classical analysis (Eq. (7)), where 
the subscript "0" means that the assumption uj <^T has been used. 

Higher order calculations of the correlation energy have been made, for example by O'Neil 



and Rostoker [^. However, in all treatments, the assumption u <^ T has been made. As 
we commented above, the assumption u; <^ T is very strong. A large part of the fluctuations 
has oj > T. 

To obtain the interaction energy, we need to subtract the energy of the particles due to 
their own fields, the second term of Eq. (5), from the longitudinal energy density, pL. We 
thus have pmt = Pl — -^ j dkh'^cp^. Using Eq. (9), the interaction energy can be written 
as Pint = |(^^)A, where pmt is the equivalent of the correlation energy. In fact, using the 
approximation u <t^T, pint is equal to the second term of Eq. (|^. 

In order to compare pint with Ec-, we define the parameter, 

f.^-^. (10) 



We previously found [Q that the transverse energy (summing the transverse electric and 
magnetic field energies, pej, and pb) has an additional energy, compared to the blackbody 
energy density in vacuum. The additional transverse energy is 

Ap-y = Pb + Pet - Pt , (11) 

where p^ is the photon energy density, estimated as the blackbody energy density in vacuum. 
Adding the interaction energy pint to Ap^, we obtain the total change in the energy 
density due to the transverse and longitudinal components, 

Pnew = ^P-y + Pint ■ (12) 

We calculate pnew and pint for an electron-proton plasma at T = 10^ K, T = 10^ K and 
T = 10^ K for densities ranging from 10^ — 10^^ cm~^. The densities were chosen so as to 
assure that the plasma parameter, g = l/n\\) < 1, in order that kinetic theory is valid. 
For these plasmas, the de Broglie wavelength is less than the distance of closest approach of 
thermal electrons, which justifies our classical treatment. 

In Figure 1, we plot A as a function of the density 10^ cm~^ ^ n < 10^^ cirT^ for the 
temperatures T = 10^ K, 10^ K, and 10^ K. We extended each plot until the density for 
which g = 0.3 was reached. For each of the temperatures, the value of g increases with the 
density. In the case of T = 10^ K, for example, for n = 10^ cm~^, g = 9.62 x 10~^ and for 
n = 10^ cm-^, g = 3.04 x 10"^. When g = 0.3, n = 10^^ cm-^. In the case of T = 10^ K, 
for n = 10^ cm-^, g = 3.04 x lO"*^ and for n = 10^° cm~^, g = 9.62 x 10"^ When g = 0.3, 
n = 10^^ cm~^. 

We found a very good fit for the results of Figure 1, using a Fermi-Dirac functional form 
for the density dependence of A, A(T) = Al / {exp[{x / A2) — A3] + 1), with x = login) and 
Al = aiQ + aiiT + a^T'^; A2 = a2o + a2iT + a22T'^ and y43 = a^o + a^iT + a32T'^ ■ From Figure 
1, we obtain Al = 0.3522 - 0.1698(T/10^) +0.1145(T/10^)^ A2 = 0.8255 + 0.4797(T/105) - 
0.4532(r/105)^ and A3 = 17.650 + 33.027(T/105) - 26.201 (T/10^)^ The curves (filled, 
dashed and dotted) are evaluated from the analytic expression; the filled squares are the 



calculated values of A from Eqs. (l)-(4). The fit can be seen to be excellent. In Figure 2, 
we plot F = {A — Ao)/Ao as a function of the density, for the temperatures T = 10^ K, 
10^ K and 10^ K, which shows how A differs from the usual correction Aq. 

The values of A that we obtained are positive and larger in absolute value than Aq, 
whereas Aq is negative. This indicates that the energy in the fluctuations dominates the 
interaction energy of the particles. We observe that F can reach values of a thousand or 
greater. 

The additional transverse energy Ap^ is completely negligible for these temperatures 
and densities. For example, for T = 10^ K and n = 10^^ cm~^, Ap = 10^'^p^. For this 
temperature and density, Ppar = 273p^ and pnew is completely dominated by pmt = Ap^ar- 
For example, for T = 10^ K and n = 10^^ cm^^, Ap^ = lO^^p^. 

As a check, we calculated pmt, integrating in frequency only up to cj = Up, the plasma 
frequency (^ T), and integrating in wavenumber up to k < ko- As expected, we then found 
that A is equal to Aq, the value obtained in previous analysis. 

IV. CONCLUSIONS AND DISCUSSION 

We calculated pnew and pint for an electron-proton plasma as a function of density for 
T = 10^ — 10^ K. For many interesting plasmas, we found that A ^ Aq. We used the BGK 
collison term as a rough guide to the inclusion of collisions. The BGK is a model collision 
term for the Boltzmann collision term. Collisions, however, change the results very little. For 
example, for T = 10^ K and n = 10^° cirT^ ^ the difference in A, with or without collisions, 
is less than 10^^. Since there is no significant difference between the energy density, with or 
without the collision term, the use of a more acurate collision term than the BGK collision 
term is not necessary. 

Appreciably different values than the usual ones are obtained, for the interaction energy 
of a plasma, by not assuming u <^T . To the first order in g, we found that the energy of 
an ideal gas needs to be corrected by a positive value, approximately O.Sppar = 0.3(3/2)nT. 

7 



This results in very different values from the usual ones ~ 10 ^ — 10 ^ (3/2)nT. 

We obtained a general expression for the correction A as a function of density and 
temperature: A(T) = Al / {exp[{x / A2) - A3] + 1) with x = log{n), Al = 0.3522 - 
0.1698(T/105) + 0.1145(T/105)^ A2 = 0.8255 + 0.4797(T/105) - 0.4532(T/105)^ and 
A3 = 17.650 + 33.027(T/105) - 26.201(T/105)^ 

The total correction to the energy is completely dominated by the interaction energy. 
For these temperatures and densities, the transverse additional energy is negligible. 

Our results may be applied to the plasma before the recombination era, when the plasma 
had a temperature T > 10'^ K and a density n > 10^ cm"^. Since the expansion rate of the 
Universe (the Hubble parameter) is proportional to the square root of the plasma energy 
density, our results indicate that the Universe before the recombination era was expanding 
appreciably faster than previously thought. 

The purpose of this work was to demonstrate that cj ^ T is an extremely strong as- 
sumption. By not making this assumption, there is a large change in the energy of the 
plasma. 
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FIGURES 
FIG. 1. The correction A as a function of density and temperature. The filled curve is for 

T = 10^ K, the dashed curve for T = 10^ K and the dotted curve for T = 10^. The curves are 

evaluated from the analytic expression and the filled squares are the calculated values from Eqs. 

(l)-(4). 

FIG. 2. The deviation of the correction, A from the usual one, Aq: F =| A | — | Aq | / | Aq | . 
The filled curve is for T = 10^ K, the dashed curve for T = lO'^ K and the dotted curve for 
T = 10^. 
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